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A computer-controlled data collection system was designed and 
implemented to analyze finite amplitude acoustic standing waves in a 
rectangular air-filled cavity with linear and wedge-shaped boundary 
perturbations. The response of a nonlinear standing wave in per- 
turbed and unperturbed cavities was studied. The experimental 
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titative comparisons with theory could not be made because the prop- 
erties of the cavity did not match those for which theoretical 
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I. INTRODUCTION 



A one-dimensional mathematical model for finite-amplitude 
standing waves in a real cavity was studied by Coppens and Sanders 
[Refs. 1, 21 in 1968. In 1974, this model was extended to three 
dimensions. Losses and resonance frequencies, inadequately 
described by the classical theory of wall losses, are obtained for the 
real cavity from measurements obtained in the linear acoustic regime. 
The research done in 1974 showed that when a cavity is driven near 
the resonance of a standing wave the nonlinear distortion primarily 
involves only the normal modes of the cavity that are “familially 
related” to the driven mode. Lane [Ref. 3], DeVall [Ref. 4], Slocum 
[Ref. 5], and Kuntsal [Ref. 6] examined the problem and showed that 
experimental results are in excellent agreement with the theory 
except when a degeneracy exists in the cavity. 

The purpose of the research reported in this thesis was to design 
and implement a computer -assisted data collection and analysis sys- 
tem for measuring the finite-amplitude distortion in a perturbed 
cavity with degeneracies. 
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II. BACKGROUND AND THEORY 



When the sound intensity in an acoustic process becomes so high 
that the principle of superposition no longer applies, we enter the 
domain of nonlinear acoustics. The limits of linear approximations in 
the equation of force, mass continuity, and state are exceeded and the 
resultant vibrational disturbance at any point in the medium will no 
longer be equal to the sum of its individual components. The essence 
of the nonlinear acoustic mechanism is that an intense wave perturbs 
the medium in which it exists and this perturbation alters the natural 
order governing the wave's behavior. The result is that segments of 
the disturbance in a warmer compression travel faster than those in a 
cooler rarefaction. When this happens, the acoustical wave alters its 
shape, steepening as it travels along. The alteration is infinitesimal for 
weak disturbances and gradually increases in importance with 
increases in both amplitude and frequency of sound. [Ref. 2] 

In 1968, Kirchoff (Ref. 7] studied the distortion of intense acous- 
tic waves in an unbounded fluid. This study was extended by Lamb 
[Ref. 81. Fay [Ref. 9], and by Keller [Ref. 10). 

Keck and Beyer [Ref. 11] approached the problem of unbounded 
propagation by using perturbation methods and Fourier series 
representations of the waveform. The model of Coppens and Sanders 
[Ref. 12] dealt with finite-amplitude standing waves in rigid walled 
cavities. They extended the Keck-Beyer perturbation approach to 
include wall losses predicted by the Rayleigh-Kirchoff theory [Ref. 7], 



and showed excellent agreement with experiments conducted in a 
rigid-walled tube at low but finite levels of excitation. The 
experimental work by Beech (Ref. 13] and by Ruff [Ref. 14] showed 
that at high excitation levels differences developed between theory 
and experiment. Experimental work by DeVall [Ref. 4] and by Kilmer 
[Ref. 15] showed that if degeneracies existed the model failed to 
account for the experimentally observed excitation of “nonfamilial” 
modes. 

The effects of various geometrical boundary perturbations on a 

rigid-walled cavity investigated by Aydin [Ref. 16] showed that the 

standing waves that exist in an ideal cavity must be corrected when 

the boundaries are irregular. The model of Coppens and Sanders is 

based on a three-dimensional, nonlinear wave equation with a 

dissipative term describing absorptive losses encountered by standing 

waves in rigid walled cavities. As stated in [Ref. 12], 

The standing waves which can be excited in any real cavity deviate 
from predictions of the linear wave equation with ideal boundary 
conditions for the following reasons: 

a The presence of boundary-layer losses at the cavity surfaces 
yields a dispersive contribution to the wave equation. 

b. Geometrical irregularities alter the effective dimensions of the 
cavity. 

Both of these mechanisms can be treated separately as long as the 
shift in frequency is so small that the actual resonances are close to 
the theoretical values resulting from the classical model. These are 
treated by assuming the dimensions are exact, and the apparent phase 
speeds are determined on that basis. 
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The resonance frequency for each standing wave is defined as 
(Ref. 6] 



the apparent phase speed appropriate for that frequency. 

The force equation, in Eulerian coordinates, can be written in the 
form 



where the operator X describes those physical processes leading to 
absorption and dispersion. If we ignore rotational effects so that 
u = V(J), restrict waves to rectangular cavities, and neglect terms of 
orders higher than M^. M(a/k), and (a/k)2, the result is a 
quadratically nonlinear wave equation [Ref. 12]. 



The standing pressure wave in a rectangular cavity has the form 
[Ref. 12] 



(On — Cn [(nx kx)^ + [ny ky)^ + (nz kz)^]^'^ 



( 1 ) 




^ + (u • V) u +^VP = ^ Xu 
ot p p 



( 2 ) 




(3) 



p = coskx X coskyy coskzZ coscot 



(4) 



13 



If the cavity is being driven near resonance, the nonresonant terms 
contribute almost nothing with respect to the resonant terms. So the 
standing wave pressure has the form of 

oo 

P = Z Pn (5) 

n=l 



where each Pn is of frequency nco and has the form [Ref. 12). 



Pn 



= MRn cosnkxX cosnkyy cosnkzZ sin(ncot + Qn) 



( 6 ) 



These pnS constitute those normal modes familially related to pi. 

The nonlinear, coupled, transcendental equations which must be 
solved for the RnS are 



Rn 



cosl 

sin] 



(Qn - 0n) = NMpQn COsQn 



j n-l 

2 Z Rj Rn-J 

j=l 



cosl 

sin] 



* (Qj - Qn-j) 



oo 



-E 

j=i 



Rn+j Rj 



cosl 

sin] 



* (9n-j - Qj) 



for all n > 1 (7) 



R = 1 for n=l 
where Qn = Arctan(-Fn) 



Fn = Q [(nco)2 - (con)2]/(nco)2 = 2Qn for « 1 (8) 



Qn 



cOn 



(9) 



Total pressure near resonance for perturbed boundaries is [Ref. 6] 
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E 



p = (n, m, 1/(0, 9) 



CoJ 



Lx 



I 



n=o 



ami An (-1)" * 



m=o 

1=0 



Snml (n,m,l/0),6 + <7nml) (10) 

where Snml - first-order perturbation correction and Onml angle 
which is defined as [Ref. 6]; 




If the cavity has perturbed walls, the solution will be in terms of a 
summation of the classical linear solution for ideal boundaries plus 
perturbation correction terms due to the irregular boundary: 

p = Po + e p' + e2 p" + ... (11) 

Since the magnitude of the boundary perturbation is kept small, sec- 
ond and higher terms in e can be considered insignificant, so that 

p = Po + € p’ (12) 

and p must satisfy the following conditions, 

□ lP = 0 (13) 
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and 



Vp • ii = 0 at X = 0, Lx [1 + ef(y,z)] (14) 

y = 0. Ly 

z = 0. Lz 

where n, the local normal to the real surface, is obtained by taking the 
gradient of the equation for the boundary [Ref. 1], 



ri = V {x - Lx [1 + ef(y.z)]} 



(15) 



Thus, to the first order in e , 



A , 9f(y.z) A , 3f(y,z) 



(16) 



and when equation (16) is used in equation (15). the result is 



^ 3f(v,z) 3p 3f(y,z) 9p 

ax " ^ ^ av av “ ^ ^ az az 



3y 3y 



dz dz 



X = Lx [1 + ef(y,z)l 



= 0 



(17) 



Applying a Taylor-series expansion for p evaluated at the real 
boundary Lx [1 + ef(y,z)] produces 



pi = p 1 

X = Lx (1 + €f(y,z)] x = Lx 



+ 



^1 

3xx = Lx 



€ Lxf(y.z) 






( 18 ) 
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Substituting equation (12) into the RHS of equation (18). taking 
the partial derivative with respect to x on both sides and keeping the 
first-order terms in e. )delds (Ref. 1] 



^1 =-^1 

= Lx U + X = Loc = Lx 



a2p, 

+ 



eLxf(y,z) + ... 



(19) 



Taking the partial derivatives with respect to y and z and using 
exactly the same procedure gives [Ref. 1] 






§£| ^ , , -r , 

= Lx [1 + efly,z)l ^x = Lx ^x = Loc ^^x = Lx 



+ lrirl €Lxf(y.z) + ... (20) 



^1 =-^1 + 6^1 

^^x = Lx (1 + efly.z)l X = Lx x = Lx 



+ 



I 

az9xx'=Lx 



6 Lxf(y.z) 



+ ... 



( 21 ) 



Substituting equations (19), (20), and (21) into equation (17) and 
keeping the first-order terms in e results in 



^1 



, r f, 3^po 3f(y.t) 9po 3f(y.z) Spo 
= U [-f(y.z) -jjj- + + -^STJ 



X - Lx 



( 22 ) 



The RHS of equation (22) can be represented as a Fourier series in 
cosines, so that p' can be expressed as a summation of normal modes. 
Hence, 



(24) 



^1 

= Lx 



OO CO 

= E 

m=0 1=0 




ami (0,m,l/w,6) 



where 



1 

^00 " LyLz J J 



(25) 



2 ^ ^ m;t , J 

= i [GIcos— ydyd, 



(26) 



2 Lu, , 

^01 = 4i;J i [Glcos^dyd, 



(27) 



3ml 



LyLz 5 



Y Y rr-i ^ 

J J [G] COS y cos z dj 



m^tO 

n^tO 



(28) 



and G is defined as 

(29) 

OX x = Lx 

In order to find the contribution to p' from each of the terms, it is 
stated [Ref. 16] that, for a cavity forced by a dynamic boundary condi- 
tion at a boundary. 

Hl P' = 0 (30) 



and the dynamic boundary condition from the (m.l)th term is 
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(31) 



^p 'ml I 

x'=L:< 



A niTc 

A cos -f — y cos 

Ly 



Lk 




0i((ot+Q) 



then 



Pml = - A 



CO ^2 



Lx 



I A 

n=0 



, ^ nr: mic Ln 

n (-1)"^ Snml COS j^x COS "i^Y cos 



. 0i(cot+Q+anmD 



where 



A if n = 0 

- [2 if n = 1,2.3.... 



(32) 



Applying this solution to equation (24) gives the complete solution 
for the first-order perturbation [Ref. 6] 



p = (n.m.l/co,6) 




&ml An (“1)*^ Snml 



(n.m.l/w.e+Onml) 



(33) 



As it is stated (Ref. 16]. if this is near a resonance, co-conml. then 
this term will dominate the summation and all other non-degenerate 
terms can be omitted. Consequently, equation (33) becomes 



p = (n.m,l/co.e) 



3.ml An (-1)^ Snml 



(n.m.l/co.e+Onml) 



(34) 
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III. EXPERIMENTAL CONSIDERATIONS 



A. COMPUTER-CONTROLLED APPARATUS 

The rectangular cavity used in this research was designed and 
constructed in 1977 by Kuntsal [Ref. 6] (Figure 1). One wall of the 
cavity could be moved and secured. The interior of the cavity was 
12.002 inches long and 2.502 inches high, while the width could be 
varied from approximately 5.50 inches to 7.00 inches, in 0.25-inch 
steps. The piston, located in a 2.25-inch diameter port in the floor of 
the cavity, was supported by an O-ring to produce a seal between the 
piston and the port. Three 1 /4-inch diameter microphone ports 
were used. A condenser microphone was mounted in a 0.24-inch 
aluminum case and an O-ring and silicon grease were used between 
this case and the ports for sealing. 

Modal designation for the lowest three modes is shown in 
Figure 2. The (020) and (100) modes are nearly degenerate. The 
(020) mode is a one-dimensional standing wave with nodes 3 and 9 
inches along the long wall of the cavity and the (100) mode is a 
standing wave with one node in the middle of the short wall. When 
the microphone is placed at position A, the (100) mode is sensed 
without contamination from the (020) mode. In position C, the (020) 
mode is sensed without contamination from the (100) mode. A block 
diagram of the apparatus is shown in Figure 3. 

The piston was driven by an M-B Electronics model EA1500 
exciter, which in turn was driven by an APS Electronic model 114 
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Figure 1 

Rectangular Cavity 
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Figure 2 

(a) Cavity Orientation 

(b) Modal Designation 




Figure 3 

Block Diagram 
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power amplifier. The driving signal was produced by an HP 3325A 
Synthesizer/Function generator. 

The piston movement was continuously analyzed by an Endevco 
model 2215 Accelerometer mounted within the piston. The output 
voltage of the microphone and output voltage of the accelerometer 
were measured by an HP 3678 multimeter. An HP 3562A Dynamic 
Signal analyzer was used to measure the harmonic distortion, A model 
5204 Lock-in Analyzer was used to measure the resonance frequen- 
cies at small signal amplitude. 

The sound pressure in the cavity was sensed by a Bruel and Kjaer 
4136 condenser microphone with matching preamplifier 2807. The 
output of this preamplifier went to three pieces of equipment: (1) an 
HP 3478 multimeter to measure the voltage level, (2) an HP 3 5 62 A 
Dynamic Signal Analyzer to measure the harmonics, and (3) a Cos 
5060A oscilloscope to monitor the undesired distortion of the accel- 
erometer signal. 

The HP 3562A Dynamic Signal Analyzer and HP 3325A Synthe- 
sizer/Function generator were controlled by an HP 87 computer. Two 
types of geometrical perturbations were used: 

a The long wall was machined at a small angle (as shown in 
Figure 4). 

b. The long wall was machined as shown in Figure 5. 

The effect of these geometrical perturbations on the resonance 
frequency of the (100) mode could be estimated by calculating an 
“effective" width Le = L - AL of the perturbed cavity where AL = V/A 
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Figure 4 

Linearly Perturbed Wall 




Figure 5 

Wedged Perturbation 

and V and A are the volume and cross-sectional area of the 
perturbation. [Ref. 1] 

B. STRENGTH PARAMETER (STRM) 

Strength of the finite amplitude interaction is characterized by 
the strength parameter STRM [Ref. 6]. 

STRM = M(3Qi where Qi is the quality factor of the driven mode 




where Sm is the microphone sensitivity. 
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The microphone output voltage vj is read on an HP 3562A 
Dynamic Signal Analyzer set at the fundamental harmonic. Micro- 
phone sensitivity was obtained with a Bruel and Kjaer model 4220 
piston phone. 

Sm = (1.56 ± 0.03) X 10-3 volt/(N/M2) 

By using p = 1.293 kg/m3 and C = 345 m/sec, we found 

STRM = 7.07 X 10-3 xVx Q. ( 35 ) 

This is the same sensitivity as is measured by Kuntsal [Ref. 6). 

C. FREQUENCY PARAMETER (FP) 

The frequency parameter is defined by 



FP = 2Qi (CO -coi)/coi (36) 

where f is a driving resonance frequency and fi is a resonance fre- 

AcOn 

quency of (010) mode. It is assumed that « 1. FP indicates 
where the driving frequency is with respect to the fundamental reso- 
nance frequency of the cavity. As IFP I increases, the nonlinear effects 
will in general become smaller for a given strength parameter. 



D. HARMONICITY COEFFICIENT (E(n)) 

The harmonicity coefficient is defined as [Ref. 6] 



E(n) = 



cOn - ncoi 

ncoi 



(37) 
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E(n) characterizes how well the modes of a given family are tuned. 

E. COMPUTER PROGRAMS 

Two computer programs were developed. The first program finds 
the resonance frequencies and calculates the quality factors (Figure 6). 
The second program drives the Dynamic Signal Analyzer HP 3562A by 
a HP 3325A Synthesizer/Function generator and finds and calculates 
harmonic content of the microphone output at each driving frequency 
(Figure 7). 
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start 





Figure 6 

* 9 ” Calculation Flow Diagram 
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Figure 7 

Power of the Harmonic Content Calculation Flow Diagram 
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IV. DATA COLLECTION PROCEDURES 



A. SYSTEM ALIGNMENT 

The primary factor in setting the lower limit to the measurement 
of harmonic content was the amount of distortion injected into the 
cavity by the piston. Misalignment of the piston within its port causes 
distortion of the piston motion which was detected by analyzing the 
harmonic content of the accelerometer output voltage. (With the pis- 
ton removed from its port, distortion levels were of the same order of 
magnitude as the general electronic noise level of the accelerometer 
system when at rest.) Before each run. the piston was aligned within 
its port so that the accelerometer output voltage contained less than 
0.1 percent second harmonic and 0.2 percent third harmonic. 

After testing various configurations, the one selected for ease of 
alignment consisted of the cavity mounted on a tripod support directly 
above the exciter, which was supported by its own base in an upright 
position. Both the cavity tripod and the exciter rested on a 5/ 16-inch 
pad of stiff acoustic rubber. Once the system was aligned, only an 
occasional minor adjustment was necessary to maintain low piston 
distortion. 

B. AUTOMATIC DATA COLLECTION 

Since the resonance frequencies of the cavity were observed to 
drift with time, the system was allowed to warm up for at least three 
hours prior to data collection. The harmonic content of the 
accelerometer output was analyzed and the piston was adjusted until 
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the second harmonic of the accelerometer output was at least 50 dB 
below that of the fundamental frequency. Figure 8 shows the drift in fj- 
is approximately linear in time. 

Data was collected using two measurement procedures in the 
computer programs: (1) a pre- and post-run infinitesimal amplitude 

measurement: and (2) the finite amplitude measurement. 

1. Pre-Run and Post-Run Infinitesimal Amplitude 
Measurements 

During the pre-run and post-run infinitesimal amplitude 
measurement, the piston was driven 80 mV (rms). Because the (100) 
and (020) modes were nearly degenerate, measurements of the prop- 
erties of the (100) mode were made with the microphone in port C 
while measurements of the (020) mode were made with the micro- 
phone in port B. The resonant frequency of the fundamental and the 
Q factor first few (0,n,0) in the modes were determined from 



, fu-fL 

Ir- 2 



(38) 




(39) 



and the harmonicity coefficients were then calculated from 




( 40 ) 
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Figure 8 

Drift in Frequency 
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where fu and fi are the frequencies of the upper and lower half-power 
points, and fn is the resonance frequency of the nth member of the 
(010) family. Figure 9 shows the precision of frequency 
measurements. 

2. The Finite-Amplitude Measurement 

All finite-amplitude measurements were made at port B. 
Throughout this measurement, Vi had to be kept constant by the 
computer program (Appendix C), which required adjustment of the 
driving voltage applied to the piston. The frequency was increased by 
the computer in 0.5-Hz steps through ±5 Hz about the resonance 
frequency of a fundamental mode. At each driving frequency, the 
harmonic content of the microphone output, usually up to the fourth, 
was measured with the HP 3562A Dynamic Signal analyzer. 

To determine the exact value of the frequency parameter, the 
value of fr at the instant of the measurement of Vn was determined 
from the first computer program. 
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Figure 9 

Precision of Frequency Measurements 
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V. RESULTS 



In this section, the theoretical predictions and experimental 
results obtained by Kuntsal [Ref. 6 ] will be qualitatively compared to 
the experimental results obtained by using the computer-controlled 
equipment. 

The information on the empirical losses and resonance frequen- 
cies is contained in the Q’s and E’s. It is important to note that, 
because V i is always kept constant as a function of frequency, the sec- 
ond harmonic distortion is maximized for co = CO2/2. The point where 
this occurs for each V2/V1 curve is indicated by an arrow with the 
label Fo 20 - At that point, the value of the frequency parameter is 

F020 = 2 Qoio E020 
Similarly, 

Fo 30 = 2 Qoio Eoso 

is indicated on each graph. The arrow labelled Fioo indicates the 
position where the nearly degenerate ( 1 , 0 , 0 ) mode is resonant, and 
the value of F 100 is 

Fioo = 2 Qoio Eioo 

The theoretical predictions for unperturbed cavities made by 
Kuntsal [Ref. 6 ] were compared to the experimental results as seen in 
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Figure A. It is clearly seen that the theory and experiment are qualita- 
tively in agreement. 

At Figure 9. fi and f 2 are plotted as a function of time. This figure 
showed two things: (1) changes of resonance frequencies of the family 
modes are proportional to the frequency, as would be expected if the 
long dimension of the cavity changes with time, most likely due to 
temperature changes; and (2) the fluctuations of resonance frequen- 
cies from measurement to measurement indicate a precision of ± 0.02 
Hz in determining the resonance frequencies of the lowest mode. By 
using the value of ± 0.02 Hz. expected uncertainties in the E’s are 

fn ± 0.02 - n(fi ± 0.02) 

^Eexpected - 

Calculated AEexpected and measured AE’s between pre- and post-run 
are tabulated in Table 4. The difference of the pre-run and post-run 
E’s for the family modes in the unperturbed cavity are within the 
expected precision. For the 100 mode (non-family mode), this differ- 
ence is much greater than expected, indicating that the cavity may be 
changing its dimensions unequally in all directions. 

For a linear perturbation (Figure B) and for the wedge perturba- 
tion (Figure C). the harmonic content in Reference 6 is in qualitative 
agreement with experiment, considering that lower Q’s and different 
E’s were observed compared to Kuntsal’s measurement [Ref. 6). 
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VI. CONCLUSIONS 



The experimental apparatus, procedures, and computer programs 
are capable of providing data sufficiently precise to test the prediction 
of the theory in the absence and presence of geometrical perturbation. 

Because the strength parameters (STRM) and quality factors (Q) 
calculated in ,the experiment did not match the STRM and Q in Refer- 
ence 6, quantitative agreement could not be tested. However, the 
experimental results show that a geometrical perturbation alters the 
finite amplitude behavior of the cavity, and that the nature of these 
changes is in qualitative agreement with the predictions of the theory. 

The quantitative disagreements between the theory and experi- 
ment for the perturbed cavities may be because the resonance fre- 
quency of the non-family mode is not tracking the family modes. 

In order to avoid disturbing the position of the piston in the port 
of the cavity, the microphone has to be shifted very carefully from one 
port to another between the “pauses” as the program goes on. 
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APPENDIX A 



CURVES 



Average Values for Qs and Es 


Mode: 


010 


Qj = 194.72 


Ej = 0.0 


STRM: 


.389 


= 278.26 


Eg = 96.08x 10-5 






Qg = 355.18 


Eg = 16.86x10-^ 






263.12 


Ejgg = 82.67 X 10-5 




Figure A 

Unperturbed Cavity 
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Average Values for Qs and Es 


Mode: 


010 


9, = 192.81 


Ej = 0.0 


STRM: 


.355 


9^ = 287.86 


= 20.39x 10-^ 






9a = 340.61 


Eg = 30.70 X 10^ 






9,00 = 262.99 


E,oo= 23.75x 10-^ 




- 3 - 2-10123 

Frequency Parameter 



Figure B 

Linear Perturbed Cavity 
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Average Values for Qs and Es 


Mode: 


010 


Q, = 187.52 


E, = 0.0 


STRM: 


.370 


= 288.13 


Ej = -49.02 X 10-^ 






Q 3 = 350.42 


E 3 = -23.25 X 10-^ 






Q.oo= 291.40 


Ej 33 = -51.34x10-^ 




Figure C 

Wedge Perturbed Cavity 
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APPENDIX B 



TABLES 



TABLE 1 

UNPERTURBED CAVITY 

(SEE FIGURE A) 



Pre-Run 



Mode 


fu 


fL 


fr 


9 


E 


010 


565.801 


562.875 


564.351 


195.40 


0.0 


020 


1131.840 


1127.791 


1129.812 


278.69 


98.16x10-5 


030 


1698.271 


1693.524 


1695.901 


356.51 


16.82x10-4 


100 


1131.872 


1127.590 


1129.731 


264.20 


91.17x10-5 



Finite-Amplitude Run 

Vi = -11 dB Mode: 010 STRM = .389 



f (Hz) 


V2 (dB) 


V3 (dB) 


V4 (dB) 


V2/V1 


V3/V1 


V4/V1 


FP 


558.5 


-50.0 


-70.5 


-70.5 


0.0113 


0.0014 


0.0012 


-4.19 


559.0 


-47.9 


-72.0 


-69.0 


0.0143 


0.0013 


0.0011 


-3.85 


559.5 


-46.9 


-73.8 


-67.8 


0.0160 


- 


0.0011 


-3.50 


560.0 


-46.0 


-72.6 


-67.5 


0.0179 


- 


0.0010 


-3.15 


560.5 


-44.8 


-71.2 


-69.7 


0.0204 


0.0011 


0.0010 


-2.81 


561.0 


-43.7 


-69.7 


-73.5 


0.0233 


0.0012 


- 


-2.47 


561.5 


-42.4 


-67.8 


-76.5 


0.0269 


0.0014 


- 


-2.12 


562.0 


-41.0 


-65.7 


-75.9 


0.0316 


0.0018 


- 


-1.78 


562.5 


-39.4 


-63.2 


-82.5 


0.0380 


0.0025 


- 


-1.43 


563.0 


-37.8 


-60.3 


-85.0 


0.0457 


0.0034 


- 


-1.09 


563.5 


-35.9 


-56.6 


-74.9 


0.0569 


0.0052 


- 


-0.74 


564.0 


-34.0 


-52.2 


-67.5 


0.0708 


0.0070 


0.0015 


-0.40 


564.5 


-32.2 


-47.6 


-60.6 


0.0876 


0.0148 


0.0033 


-0.05 


565.0 


-31.1 


-45.0 


-54.5 


0.0994 


0.0200 


0.0067 


0.29 


565.5 


-31.6 


-43.9 


-52.9 


0.0933 


0.0226 


0.0070 


0.64 


566.0 


-33.3 


-46.7 


-56.0 


0.0772 


0.0164 


0.0056 


0.98 


566.5 


-35.3 


-51.0 


-63.2 


0.0610 


0.0100 


0.0025 


1.33 


567.0 


-37.6 


-55.5 


-70.8 


0.0484 


0.0060 


0.0010 


1.67 


567.5 


-39.8 


-60.0 


-73.6 


0.0365 


0.0036 


- 


2.02 


568.0 


-42.1 


-64.3 


-74.7 


0.0279 


0.0022 


- 


2.36 


568.5 


-44.1 


-68.5 


-74.2 


0.0221 


0.0013 


- 


2.71 


569.0 


-46.2 


-69.5 


- 


0.0174 


0.001 1 


- 


3.05 
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o 



6 



12 " 



Post-Run 



Mode 


fu 


fL 


fr 


9 


E 


010 


566.381 


563.270 


564.82 


194.03 


0.0 


020 


1132.740 


1128.676 


1130.71 


277.82 


94.00x10-5 


030 


1699.770 


1694.974 


1697.37 


353.84 


16.90x10-4 


100 


1132.800 


1128.487 


1130.64 


262.03 


74.16x10-5 
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TABLE 2 



LINEAR PERTURBATION CONFIGURATION 

(SEE FIGURE B) 



Pre-Run 



Mode 


fu 


fL 


fr 


9 


E 


010 

020 

030 

100 


569.465 

1133.883 

1703.014 

1131.650 


566.545 

1129.961 

1698.008 

1131.341 


565.653 

1132.306 

1700.594 

1134.651 


193.72 

288.71 

339.72 
263.321 


0.0 

17.67x10-4 

21.42x10-4 

29.57x10-4 



Finite-Amplitude Run 

Vi = -12 dB Mode: 010 STRM = .355 



f (Hz) 


V 2 (dB) 


V 3 (dB) 


V 4 (dB) 


V 2 /V 1 


V 3 /V 1 


V 4 /V 1 


FP 


561.0 


-43.0 


-68.7 




0.0251 


0.0013 




-3.19 


561.5 


-42.7 


- 68.0 


- 


0.0260 


0.0014 


- 


-2.84 


562.0 


-42.3 


-66.9 


- 


0.0271 


0.0016 


- 


-2.50 


562.5 


-41.6 


-64.5 


- 


0.0296 


0.0021 


- 


-2.16 


563.0 


-40.3 


-62.7 


- 


0.0342 


0.0026 


- 


-1.82 


563.5 


-38.7 


-60.1 


- 


0.0411 


0.0035 


- 


-1.47 


564.0 


-37.3 


-56.2 


- 


0.0485 


0.0055 


- 


-1.13 


564.5 


-34.2 


-53.5 


- 


0.0691 


0.0075 


- 


-0.79 


565.0 


-33.8 


-51.0 


-65.0 


0.0721 


0.0100 


0.0020 


-0.44 


565.5 


-31.8 


-46.4 


-60.1 


0.0911 


0.0170 


0.0035 


- 0.10 


566.0 


-30.9 


-43.0 


-54.6 


0.1012 


0.0251 


0.0066 


0.24 


566.5 


-29.6 


-40.3 


-50.0 


0.1161 


0.0344 


0.0112 


0.58 


567.0 


-30.9 


-45.8 


-54.1 


0.1000 


0.0180 


0.0070 


0.92 


567.5 


-32.8 


-49.8 


-59.2 


0.0810 


0.0115 


0.0039 


1.26 


568.0 


-34.1 


-51.5 


-65.0 


0.0700 


0.0094 


0.0020 


1.61 


568.5 


-36.7 


-56.2 


- 68.1 


0.0514 


0.0055 


0.0014 


1.95 


569.0 


-38.0 


-58.7 


- 


0.0441 


0.0041 


- 


2.29 
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Post-Run 



Mode 


fu 


fL 


fr 


9 


E 


010 


570.410 


567.456 


566.861 


191.90 


0.0 


020 


1135.067 


1131.116 


1133.984 


287.01 


23.10x10-4 


030 


1705.321 


1700.340 


1701.110 


341.50 


39.98x10-4 


100 


1136.740 


1132.416 


1135.756 


262.66 


17.94x10-4 
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TABLE 3 



WEDGE PERTURBATION CONFIGURATION 

(SEE FIGURE C) 



Pre-Run 



Mode 


fu 


fL 


fr 


Q 


E 


010 

020 

030 

100 


566.912 

1127.191 

1695.310 

1126.761 


563.864 

1123.281 

1690.471 

1122.911 


565.381 

1125.230 

1692.891 

1124.829 


185.80 

288.08 

349.27 

291.41 


0.0 

-48.92x10-4 

-19.17x10-4 

-52.47x10-4 



Finite-Amplitude Run 

Vi = -11 dB Mode: 010 STRM = .370 



f (Hz) 


V2 (dB) 


V3 (dB) 


V4 (dB) 


V2/V1 




V4/V1 


FP 


559.4 


-39.5 


-66.9 


_ 


0.0376 


0.0016 


_ 


-3.99 


559.9 


-37.9 


-65.3 


- 


0.0450 


0.0019 


- 


-3.66 


560.4 


-36.4 


-63.9 


- 


0.0540 


0.0023 


- 


-3.33 


560.9 


-34.5 


-61.8 


- 


0.0672 


0.0029 


- 


-3.00 


561.4 


-32.4 


-59.1 


-78.0 


0.0856 


0.0039 


- 


-2.66 


561.9 


-30.3 


-55.7 


-72.7 


0.1090 


0.0058 


0.0008 


-2.33 


562.4 


-28.7 


-55.6 


-68.2 


0.1310 


0.0059 


0.0014 


-2.00 


562.9 


-28.9 


-50.8 


-65.3 


0.1271 


0.0103 


0.0017 


-1.67 


563.4 


-30.9 


-50.0 


-62.5 


0.1010 


0.0112 


0.0027 


-1.34 


563.9 


-32.3 


-49.7 


-60.3 


0.0810 


0.0116 


0.0034 


-1.01 


564.4 


-34.7 


-50.9 


-61.4 


0.0650 


0.0101 


0.0030 


-0.67 


564.9 


-36.3 


-54.3 


-67.0 


0.0550 


0.0068 


0.0016 


-0.34 


565.4 


-37.5 


-59.3 


-75.2 


0.0480 


0.0038 


- 


-0.01 


565.9 


-38.4 


-65.0 


- 


0.0430 


0.0020 


- 


0.32 


566.4 


-39.3 


-69.7 


- 


0.0391 


0.0012 


- 


0.65 


566.9 


-39.7 


-72.2 


- 


0.0370 


0.0009 


- 


0.98 


567.4 


-40.6 


-74.5 


- 


0.0280 


- 


- 


1.32 


567.9 


-40.7 


-75.4 


- 


0.0230 


- 


- 


1.65 


568.4 


-41.6 


-73.5 


- 


0.0210 


- 


- 


1.98 


568.9 


-44.0 


-71.5 


- 


0.0193 


- 


- 


2.31 


569.4 


-47.3 


-70.9 


- 


0.0151 


- 


- 


2.64 


569.9 


-49.1 


-71.1 


- 


0.0120 


- 


- 


2.97 
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Post-Run 



Mode 


fu 


fL 


fr 


9 


E 


010 


566.950 


563.963 


565.464 


189.24 


0.0 


020 


1127.311 


1123.411 


1125.371 


288.18 


-49.12x10-4 


030 


1695.544 


1690.731 


1693.144 


351.57 


-27.33x10-4 


100 


1 126.981 


1123.126 


1125.051 


291.39 


-50.21x10-4 
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TABLE 4 



E’S DIFFERENCES BETWEEN RUNS 





Unperturbed 


Linear Perturbed 


Wedge Perturbed 




Configuration 


Configuration 


Configuration 


AE2 


.416x10-4 


5.43x10-4 


0.2x10-4 


AE 3 


.08x10-4 


8.56x10-4 


8.16x10-4 


AEioo 


1.7x10-4 


11.63x10-4 


2.26x10-4 


AEexpected 


.53x10-4 
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APPENDIX C 



COMPUTER PROGRAMS 



This program drives the HP 3325A Synthesizer/ Function genera- 
tor between 562 Hz and 1710 Hz. finds the resonance frequency, and 
calculates the quality factor (Q). 



10 ! SUBROUTINE Q 

14 SETTIME 0,87200 

1 5 FOR INC=1 TO 30 
20 K=-1 

30 FORFREQ=562TO1710STEP2 

50 IF FREO>550 AND FREQ<570 THEN GOTO 21 0 

70 IF FREO>570 THEN GOTO 90 

90 FOR FREO=1116TO 1140STEP2 

110 IF FREO>1 1 1 5 AND FREQ<1 1 40 THEN GOTO 21 0 

130 IF FRO>1 140 THEN GOTO 1791 

140 K=K+1 

150 FORFREO=1680TO1710STEP2 

1 70 IF FREO>1 670 AND FREQ<1 71 0 THEN GOTO 21 0 

190 IF FREO>1 710 THEN GOTO 1800 

210 F1=FREO 

230 F2=F1+2 

250 F3=F1+4 

270 REMOTE 717,723,702 

290 OUTPUT 71 7 ;"FR";F1 ;"HZ" 

310 WAIT 5000 

330 ENTER 723 ;AMP1 

350 OUTPUT 71 7 ;"FR";F2;"HZ 

370 WAIT 5000 

390 ENTER 723 ;AMP2 

41 0 OUTPUT 71 7 ;"FR";F3;"HZ" 

430 WAIT 5000 

450 ENTER 723 ;AMP3 

470 DISP FI ,F2,F3,AMP1 ,AMP2,AMP3,K 

51 0 IF AMP2<AMP1 OR AMP2<AMP3 THEN GOTO 590 

570 IF AMP2>AMP1 OR AMP2>AMP3 THEN GOSUB BINOM 

590 NEXT FREQ 

610 ! SUBROUTINE BINOM 
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630 BINOM: 

650 P=F2 
670 S=2 

690 S=S/2 

71 0 OUTPUT 71 7 ;"FR";P;"HZ" 

720 WAIT 5000 

730 ENTER 723 ; AP 

750 IF P>1 710 THEN GOTO 1800 

770 PU=P+S 

790 OUTPUT 71 7 ;"FR";PU:"HZ" 

800 WAIT 5000 
810 ENTER 723 ;AU 
830 PL=P-S 

850 OUTPUT 71 7 ;"FR";PL;"HZ" 

860 WAIT 5000 
870 ENTER 723 ; AL 
890 DISP P.PU,PL,AP,AU,AL 
910 IF AU>APTHEN P=PU 
930 IF AL>AP THEN P=PL 
950 IF S>.008 THEN GOTO 690 
970 GOSUB HPWR 
990 RETURN 

1010 ! SUBROUTINE HALFPOWER 
1030 HPWR; 

1050 DISPP.AP 
1060 DISP TIME 
1065 PRINT TIME 
1070 PRINT "MAX" 

1090 PRINT P,AP 

1100 IF INC>1 AND INC<30 THEN GOTO 1 771 ELSE GOTO 1110 
1110 PRINT "REAL HALFPOWER" 

1130 APHl=AP/2'^.5 
1150 PRINT APH1 
1 1 70 FOR F=1 TO 2 
1190 PH=P+((-1)^F) 

1210 OUTPUT 71 7 ;"FR";PH;"HZ" 

1230 WAIT 5000 
1250 ENTER 723; AH 
1270 DISP PH, AH 

1290 IFAH<APH1 THEN GOTO 1350 ELSE GOTO 1310 
1310 PH=PH+((-1)''F) 

1330 GOTO 1210 
1350 S=.5 

1370 PHI=PH-(S*((-1)^F)) 

1390 FOR 1=1 TO 10 
1 41 0 OUTPUT 71 7 ;"FR";PHI;"HZ" 

1430 WAIT 5000 
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1450 ENTER 723 ;AHI 
1470 S=S/2 

1 490 IF AHI>APH1 THEN GOTO 1570 
1510 PHI=PHI-(S*((-1)^F)) 

1530 NEXT I 

1550 GOTO 1610 

1570 PHI=PHI+(S*((-1)''F)) 

1590 NEXT I 

1610 PRINT "HALFPOWER" 

1630 DISPPHLAHI 

1650 PRINT PHI, AHI 

1670 IF F=1 THEN PHIR=PHI 

1960 NEXTF 

1710 0=P/(PHI-PHIR) 

1730 PRINT "O VALUE" 

1750 PRINT Q 

1770 DISPO 

1771 K=K+1 

1772 IF K=0 THEN GOTO 90 

1 773 ! AFTER FIRST PAUSE SHIFT THE MICROPHONE TO THE PORT C TO 
MEASURE THE FREQUENCY OF THE DEGENERATE MODE. 

1 775 ! AFTER SECOND PAUSE SHIFT THE MICROPHONE TO THE PORT B 

1791 IF K=1 THEN PAUSE 

1792 IF K=2 THEN PAUSE 

1 793 IF K=1 THEN GOTO 90 

1794 IF K=2 THEN GOTO 140 
1797 RETURN 

1800 NEXT INC 
1810 END 
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This program drives the HP 3562A Dynamic SignalAnalyzer by HP 
3325A Synthesizer/Function generator and finds and computes the 
harmonic content of the microphone output at each driving frequency 
between 563 Hz and 570 Hz. 



10 REMOTE 720,717 

1 5 SETTIME 0,87033 

20 OUTPUT 720 :"PRST" 

30 OUTPUT 720 ;"B" 

40 OUTPUT 720 :"FRS3KHZ" 

50 OUTPUT 720 ;"PSPC' 

60 P=800 

70 OUTPUT 71 7 ;"AMP";P;"MV" 

80 OUTPUT 720 ;"LNRS" 

1 00 OUTPUT 720 ;”STBL" 

1 1 0 OUTPUT 720 ;"PSP2" 

1 20 OUTPUT 720 :"KEY31 " 

1 25 ! DRIVING THE DYNAMIC SIGNAL ANALYZER AND READING THE 

HARMONIC CONTENT OF THE DRIVING FREQUENCY. 

1 30 FREQ-563 

1 40 OUTPUT 71 7 :"FR";FREQ;"HZ" 

1 50 OUTPUT 720 SF563HZ " 

160 WAIT 3000 
1 70 OUTPUT 720 ;"STRT" 

1 80 WAIT 8000 
1 90 OUTPUT 720 ;"MKPK" 

200 OUTPUT 720 ;"RDMK" 

210 ENTER 720 : XI, YB 
220 DIL=YB 
230 DISPX1,DIL 

234 PRINT TIME 

235 PRINT XI ,DIL 

240 FOR FREQ=563 TO 570 STEP .5 

250 OUTPUT 720 :"FRS3KHZ" 

251 OUTPUT 720 :"SF";FREQ;"HZ" 

255 OUTPUT 71 7 ;"FR";FREQ:"HZ" 

260 FOR l=1 TO 4 

270 OUTPUT 71 7 :"AMP":P;"MV" 

280 OUTPUT 720 r'STRT’ 

290 WAIT 8000 

300 OUTPUT 720 ;"MKPK" 
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310 

320 

324 

325 

326 

332 

333 

335 

350 

351 

355 

356 

360 

365 

370 

380 

385 

390 

400 

410 

420 

430 

440 

450 

460 

470 

475 

480 

490 

510 

660 



OUTPUT 720 :"RDMK" 

ENTER 720 ; X.YB 
DISP X.YB 
PRINT X.YB 
PRINT TIME 

IF 1=1 THEN GOTO 333 ELSE GOTO 350 
IF FREQ>563 THEN GOTO 335 ELSE GOTO 350 
GOSUB BINOM 
OUTPUT 720 :"SF":X+50:"HZ" 

OUTPUT 720 FRS.7KHZ" 

NEXT I 
NEXT FREQ 
! SUBROUTINE BINOM 

! THIS SUBROUTINE KEEPS THE FUNDAMENTAL CONSTANT. 
BINOM: 

RAT=1 0^(DIU1 0)/1 0^(YB/1 0) 

DISP RAT.DIL.YB 

IF RAT>.999 AND RAT<1.001 THEN RETURN 
P=P*RAT\2 

OUTPUT 717 ;"AMP";P;"MV" 

OUTPUT 720 :"STRT" 

WAIT 8000 

OUTPUT 720 :"MKPK" 

OUTPUT 720 ;"RDMK" 

ENTER 720 ; X.YB 
DISP X.YB 
PRINT TIME 
PRINT X.YB 

IF FREQ>570 GOTO 660 

GOTO 380 

END 
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